hoxc12/c13 as key regulators for rebooting the developmental program in Xenopus limb regeneration.

During organ regeneration, after the initial responses to injury, gene expression patterns similar to those in normal development are reestablished during subsequent morphogenesis phases. This supports the idea that regeneration recapitulates development and predicts the existence of genes that reboot the developmental program after the initial responses. However, such rebooting mechanisms are largely unknown. Here, we explore core rebooting factors that operate during Xenopus limb regeneration. Transcriptomic analysis of larval limb blastema reveals that hoxc12/c13 show the highest regeneration specificity in expression. Knocking out each of them through genome editing inhibits cell proliferation and expression of a group of genes that are essential for development, resulting in autopod regeneration failure, while limb development and initial blastema formation are not affected. Furthermore, the induction of hoxc12/c13 expression partially restores froglet regenerative capacity which is normally very limited compared to larval regeneration. Thus, we demonstrate the existence of genes that have a profound impact alone on rebooting of the developmental program in a regeneration-specific manner.

Fgf10 mutant newts regenerate normal hindlimbs despite severe developmental defects.

In amniote limbs, Fibroblast Growth Factor 10 (FGF10) is essential for limb development, but whether this function is broadly conserved in tetrapods and/or involved in adult limb regeneration remains unknown. To tackle this question, we established Fgf10 mutant lines in the newt Pleurodeles waltl which has amazing regenerative ability. While Fgf10 mutant forelimbs develop normally, the hindlimbs fail to develop and downregulate FGF target genes. Despite these developmental defects, Fgf10 mutants were able to regenerate normal hindlimbs rather than recapitulating the embryonic phenotype. Together, our results demonstrate an important role for FGF10 in hindlimb formation, but little or no function in regeneration, suggesting that different mechanisms operate during limb regeneration versus development.

A CRISPR-Cas9-mediated versatile method for targeted integration of a fluorescent protein gene to visualize endogenous gene expression in Xenopus laevis.

Xenopus laevis is a widely used model organism in developmental and regeneration studies. Despite several reports regarding targeted integration techniques in Xenopus, there is still room for improvement of them, especially in creating reporter lines that rely on endogenous regulatory enhancers/promoters. We developed a CRISPR-Cas9-based simple method to efficiently introduce a fluorescent protein gene into 5' untranslated regions (5'UTRs) of target genes in Xenopus laevis. A donor plasmid DNA encoding an enhanced green fluorescent protein (eGFP) flanked by a genomic fragment ranging from 66 bp to 878 bp including target 5'UTR was co-injected into fertilized eggs with a single guide RNA and Cas9 protein. Injections for krt12.2.L, myod1.S, sox2.L or brevican.S resulted in embryos expressing eGFP fluorescence in a tissue-specific manner, recapitulating endogenous expression of target genes. Integrations of the donor DNA into the target regions were examined by genotyping PCR for the eGFP-expressing embryos. The rate of embryos expressing the specific eGFP varied from 2.1% to 13.2% depending on the target locus and length of the genomic fragment in the donor plasmids. Germline transmission of an integrated DNA was observed. This simple method provides a powerful tool for exploring gene expression and function in developmental and regeneration research in X. laevis.

Protocols for transgenesis at a safe harbor site in the Xenopus laevis genome using CRISPR-Cas9.

We have established a new transgenesis protocol based on CRISPR-Cas9, "New and Easy XenopusTransgenesis (NEXTrans)," and identified a novel safe harbor site in African clawed frogs, Xenopus laevis. We describe steps in detail for the construction of NEXTrans plasmid and guide RNA, CRISPR-Cas9-mediated NEXTrans plasmid integration into the locus, and its validation by genomic PCR. This improved strategy allows us to simply generate transgenic animals that stably express the transgene. For complete details on the use and execution of this protocol, please refer to Shibata et al. (2022).1.

CRISPR-Cas9-Based Functional Analysis in Amphibians: Xenopus laevis, Xenopus tropicalis, and Pleurodeles waltl.

Amphibians have made many fundamental contributions to our knowledge, from basic biology to biomedical research on human diseases. Current genome editing tools based on the CRISPR-Cas system enable us to perform gene functional analysis in vivo, even in non-model organisms. We introduce here a highly efficient and easy protocol for gene knockout, which can be used in three different amphibians seamlessly: Xenopus laevis, Xenopus tropicalis, and Pleurodeles waltl. As it utilizes Cas9 ribonucleoprotein complex (RNP) for injection, this cloning-free method enables researchers to obtain founder embryos with a nearly complete knockout phenotype within a week. To evaluate somatic mutation rate and its correlation to the phenotype of a Cas9 RNP-injected embryo (crispant), we also present accurate and cost-effective genotyping methods using pooled amplicon-sequencing and a user-friendly web-based tool.

CRISPR/Cas9-based simple transgenesis in Xenopus laevis.

Transgenic techniques have greatly increased our understanding of the transcriptional regulation of target genes through live reporter imaging, as well as the spatiotemporal function of a gene using loss- and gain-of-function constructs. In Xenopus species, two well-established transgenic methods, restriction enzyme-mediated integration and I-SceI meganuclease-mediated transgenesis, have been used to generate transgenic animals. However, donor plasmids are randomly integrated into the Xenopus genome in both methods. Here, we established a new and simple targeted transgenesis technique based on CRISPR/Cas9 in Xenopus laevis. In this method, Cas9 ribonucleoprotein (RNP) targeting a putative harbor site (the transforming growth factor beta receptor 2-like (tgfbr2l) locus) and a preset donor plasmid DNA were co-injected into the one-cell stage embryos of X. laevis. Approximately 10% of faithful reporter expression was detected in F0 crispants in a promoter/enhancer-specific manner. Importantly, efficient germline transmission and stable transgene expression were observed in the F1 offspring. The simplicity of this method only required preparation of a donor vector containing the tgfbr2l genome fragment and Cas9 RNP targeting this site, which are common experimental procedures used in Xenopus laboratories. Our improved technique allows the simple generation of transgenic X. laevis, so is expected to become a powerful tool for reporter assay and gene function analysis.

Newt Hoxa13 has an essential and predominant role in digit formation during development and regeneration.

The 5'Hox genes play crucial roles in limb development and specify regions in the proximal-distal axis of limbs. However, there is no direct genetic evidence that Hox genes are essential for limb development in non-mammalian tetrapods or for limb regeneration. Here, we produced single to quadruple Hox13 paralog mutants using the CRISPR/Cas9 system in newts (Pleurodeles waltl), which have strong regenerative capacities, and also produced germline mutants. We show that Hox13 genes are essential for digit formation in development, as in mice. In addition, Hoxa13 has a predominant role in digit formation, unlike in mice. The predominance is probably due to the restricted expression pattern of Hoxd13 in limb buds and the strong dependence of Hoxd13 expression on Hoxa13. Finally, we demonstrate that Hox13 genes are also necessary for digit formation in limb regeneration. Our findings reveal that the general function of Hox13 genes is conserved between limb development and regeneration, and across taxa. The predominance of Hoxa13 function both in newt limbs and fish fins, but not in mouse limbs, suggests a potential contribution of Hoxa13 function in fin-to-limb transition.

Fgf10-CRISPR mosaic mutants demonstrate the gene dose-related loss of the accessory lobe and decrease in the number of alveolar type 2 epithelial cells in mouse lung.

CRISPR/Cas9-mediated gene editing often generates founder generation (F0) mice that exhibit somatic mosaicism in the targeted gene(s). It has been known that Fibroblast growth factor 10 (Fgf10)-null mice exhibit limbless and lungless phenotypes, while intermediate limb phenotypes (variable defective limbs) are observed in the Fgf10-CRISPR F0 mice. However, how the lung phenotype in the Fgf10-mosaic mutants is related to the limb phenotype and genotype has not been investigated. In this study, we examined variable lung phenotypes in the Fgf10-targeted F0 mice to determine if the lung phenotype was correlated with percentage of functional Fgf10 genotypes. Firstly, according to a previous report, Fgf10-CRISPR F0 embryos on embryonic day 16.5 (E16.5) were classified into three types: type I, no limb; type II, limb defect; and type III, normal limbs. Cartilage and bone staining showed that limb truncations were observed in the girdle, (type I), stylopodial, or zeugopodial region (type II). Deep sequencing of the Fgf10-mutant genomes revealed that the mean proportion of codons that encode putative functional FGF10 was 8.3 ± 6.2% in type I, 25.3 ± 2.7% in type II, and 54.3 ± 9.5% in type III (mean ± standard error of the mean) mutants at E16.5. Histological studies showed that almost all lung lobes were absent in type I embryos. The accessory lung lobe was often absent in type II embryos with other lobes dysplastic. All lung lobes formed in type III embryos. The number of terminal tubules was significantly lower in type I and II embryos, but unchanged in type III embryos. To identify alveolar type 2 epithelial (AECII) cells, known to be reduced in the Fgf10-heterozygous mutant, immunostaining using anti-surfactant protein C (SPC) antibody was performed: In the E18.5 lungs, the number of AECII was correlated to the percentage of functional Fgf10 genotypes. These data suggest the Fgf10 gene dose-related loss of the accessory lobe and decrease in the number of alveolar type 2 epithelial cells in mouse lung. Since dysfunction of AECII cells has been implicated in the pathogenesis of parenchymal lung diseases, the Fgf10-CRISPR F0 mouse would present an ideal experimental system to explore it.

Amiodarone bioconcentration and suppression of metamorphosis in Xenopus.

Trace concentrations of a number of pharmaceutically active compounds have been detected in the aquatic environment in many countries, where they are thought to have the potential to exert adverse effects on non-target organisms. Amiodarone (AMD) is one such high-risk compound commonly used in general hospitals. AMD is known to alter normal thyroid hormone (TH) function, although little information is available regarding the specific mechanism by which this disruption occurs. Anuran tadpole metamorphosis is a TH-controlled developmental process and has proven to be useful as a screening tool for environmental pollutants suspected of disrupting TH functions. In the present study, our objective was to clarify the effects of AMD on Xenopus metamorphosis as well as to assess the bioconcentration of this pharmaceutical in the liver. We found that AMD suppressed spontaneous metamorphosis, including tail regression and hindlimb elongation in pro-metamorphic stage tadpoles, which is controlled by endogenous circulating TH, indicating that AMD is a TH antagonist. In transgenic X. laevis tadpoles carrying plasmid DNA containing TH-responsive element (TRE) and a 5'-upstream promoter region of the TH receptor (TR) ßA1 gene linked to a green fluorescent protein (EGFP) gene, triiodothyronine (T3) exposure induced a strong EGFP expression in the hind limbs, whereas the addition of AMD to T3 suppressed EGFP expression, suggesting that this drug interferes with the binding of T3 to TR, leading to the inhibition of TR-mediated gene expression. We also found AMD to be highly bioconcentrated in the liver of pro-metamorphic X. tropicalis tadpoles, and we monitored hepatic accumulation of this drug using mass spectrometry imaging (MSI). Our findings suggest that AMD imposes potential risk to aquatic wildlife by disrupting TH homeostasis, with further possibility of accumulating in organisms higher up in the food chain.

A simple and practical workflow for genotyping of CRISPR-Cas9-based knockout phenotypes using multiplexed amplicon sequencing.

Founder animals carrying high proportions of somatic mutation induced by CRISPR-Cas9 enable a rapid and scalable strategy for the functional screening of numerous target genes in vivo. In this functional screening, genotyping using pooled amplicons with next-generation sequencing is the most suitable approach for large-scale management of multiple samples and accurate evaluation of the efficiency of Cas9-induced somatic mutations at target sites. Here, we present a simple workflow for genotyping of multiple CRISPR-Cas9-based knockout founders by pooled amplicon sequencing. Using custom barcoded primers, pooled amplicons from multiple individuals can be run in a single-indexed library on the Illumina MiSeq platform. Additionally, a user-friendly web tool, CLiCKAR, is available to simultaneously perform demultiplexing of pooled sequence data and evaluation of somatic mutation in each phenotype. CLiCKAR provides users with practical reports regarding the positions of insertions/deletions, as well as the frameshift ratio and tables containing mutation sequences, and read counts of each phenotype, with just a few clicks by the implementation of demultiplexing for pooled sample data and calculation of the frameshift ratio. This genotyping workflow can be harnessed to evaluate genotype-phenotype correlations in CRISPR-Cas9-based loss-of-function screening of numerous target genes in various organisms.

In vitro assessment of effects of persistent organic pollutants on the transactivation of estrogen receptor α and ß (ERα and ERß) from the Baikal seal (Pusa sibirica).

To assess the effect of exposure to persistent organic pollutants (POPs) on the estrogen receptor (ER) signaling pathway in Baikal seals (Pusa sibirica), we investigated the molecular characterizations and functions of two Baikal seal ER (bsER) isoforms, bsERα and bsERß. The bsERα and bsERß cDNA clones isolated have an open reading frame of 595 and 530 amino acid residues, respectively. The tissue distribution analyses of bsER mRNAs showed that bsERα transcripts were primarily found in the ovary and uterus, and bsERß in the muscle in wild Baikal seals. The immunofluorescence staining assay showed that 17ß-estradiol (E2) treatment promoted the nuclear translocation of in vitro-expressed bsERα. Transient transfection of bsERα in U2OS cells enhanced the transcription of pS2, an ER target gene of E2. We then measured bsER-mediated transactivation potencies of POPs in an in vitro reporter gene assay system, in which a bsERα or bsERß expression vector was transfected into COS-1 cells. For comparison, transactivation potencies of POPs on mouse ERs (mERα and mERß) were also evaluated in the same manner. Results showed significant dose-dependent responses of bsERs and mERs when treated with p,p'-dichlorodiphenyltrichloroethane (p,p'-DDT), and p,p'-dichlorodiphenyldichloroethylene (p,p'-DDE). bsERs and mERs showed no response when exposed to polychlorinated biphenyls (PCBs) or 2,3,7,8-tetrachlorodibenzo-p-dioxin. Comparison of the dose-response curves of DDTs across species (bsERs vs. mERs) showed that bsERα had a response similar to mERα, but bsERß was less sensitive than mERß. Comparing the lowest observable effective concentrations of p,p'-DDT (2.8 µM) and p,p'-DDE (10 µM) for in vitro bsERα-mediated transactivation with their hepatic concentrations in wild Baikal seals indicated that some individuals accumulated these compounds at levels comparable to the effective concentrations, suggesting the potential disruption of the bsERα signaling pathway in the wild population by these compounds. Co-transfection experiments with bsER and the aryl hydrocarbon receptor (AHR) suggested that high accumulation of estrogenic compounds exerts a synergistic effect with dioxin-like congeners on ER signaling through AHR activation in the wild seal population.

A comprehensive reference transcriptome resource for the Iberian ribbed newt Pleurodeles waltl, an emerging model for developmental and regeneration biology.

Urodele newts have unique biological properties, notably including prominent regeneration ability. The Iberian ribbed newt, Pleurodeles waltl, is a promising model amphibian distinguished by ease of breeding and efficient transgenic and genome editing methods. However, limited genetic information is available for P. waltl. We conducted an intensive transcriptome analysis of P. waltl using RNA-sequencing to build and annotate gene models. We generated 1.2 billion Illumina reads from a wide variety of samples across 12 different tissues/organs, unfertilized egg, and embryos at eight different developmental stages. These reads were assembled into 1,395,387 contigs, from which 202,788 non-redundant ORF models were constructed. The set is expected to cover a large fraction of P. waltl protein-coding genes, as confirmed by BUSCO analysis, where 98% of universal single-copy orthologs were identified. Ortholog analyses revealed the gene repertoire evolution of urodele amphibians. Using the gene set as a reference, gene network analysis identified regeneration-, developmental-stage-, and tissue-specific co-expressed gene modules. Our transcriptome resource is expected to enhance future research employing this emerging model animal for regeneration research as well as for investigations in other areas including developmental biology, stem cell biology, and cancer research. These data are available via our portal website, iNewt (http://www.nibb.ac.jp/imori/main/).

Cas9 ribonucleoprotein complex allows direct and rapid analysis of coding and noncoding regions of target genes in Pleurodeles waltl development and regeneration.

Newts have remarkable ability to regenerate their organs and have been used in research for centuries. However, the laborious work of breeding has hampered reverse genetics strategies in newt. Here, we present simple and efficient gene knockout using Cas9 ribonucleoprotein complex (RNP) in Pleurodeles waltl, a species suitable for regenerative biology studies using reverse genetics. Most of the founders exhibited severe phenotypes against each target gene (tyrosinase, pax6, tbx5); notably, all tyrosinase Cas9 RNP-injected embryos showed complete albinism. Moreover, amplicon sequencing analysis of Cas9 RNP-injected embryos revealed virtually complete biallelic disruption at target loci in founders, allowing direct phenotype analysis in the F0 generation. In addition, we demonstrated the generation of tyrosinase null F1 offspring within a year. Finally, we expanded this approach to the analysis of noncoding regulatory elements by targeting limb-specific enhancer of sonic hedgehog, known as the zone of polarizing activity regulatory sequence (ZRS; also called MFCS1). Disruption of ZRS led to digit deformation in limb regeneration. From these results, we are confident that this highly efficient gene knockout method will accelerate gene functional analysis in the post-genome era of salamanders.

Functional analysis of thyroid hormone receptor beta in Xenopus tropicalis founders using CRISPR-Cas.

Amphibians provide an ideal model to study the actions of thyroid hormone (TH) in animal development because TH signaling via two TH receptors, TRα and TRß, is indispensable for amphibian metamorphosis. However, specific roles for the TRß isoform in metamorphosis are poorly understood. To address this issue, we generated trß-disrupted Xenopus tropicalis tadpoles using the CRISPR-Cas system. We first established a highly efficient and rapid workflow for gene disruption in the founder generation (F0) by injecting sgRNA and Cas9 ribonucleoprotein. Most embryos showed severe mutant phenotypes carrying high somatic mutation rates. Utilizing this founder analysis system, we examined the role of trß in metamorphosis. trß-disrupted pre-metamorphic tadpoles exhibited mixed responsiveness to exogenous TH. Specifically, gill resorption and activation of several TH-response genes, including trß itself and two protease genes, were impaired. However, hind limb outgrowth and induction of the TH-response genes, klf9 and fra-2, were not affected by loss of trß Surprisingly, trß-disrupted tadpoles were able to undergo spontaneous metamorphosis normally, except for a slight delay in tail resorption. These results indicate TRß is not required but contributes to the timing of resorptive events of metamorphosis.

Dynamic changes in the interchromosomal interaction of early histone gene loci during development of sea urchin.

The nuclear positioning and chromatin dynamics of eukaryotic genes are closely related to the regulation of gene expression, but they have not been well examined during early development, which is accompanied by rapid cell cycle progression and dynamic changes in nuclear organization, such as nuclear size and chromatin constitution. In this study, we focused on the early development of the sea urchin Hemicentrotus pulcherrimus and performed three-dimensional fluorescence in situ hybridization of gene loci encoding early histones (one of the types of histone in sea urchin). There are two non-allelic early histone gene loci per sea urchin genome. We found that during the morula stage, when the early histone gene expression levels are at their maximum, interchromosomal interactions were often formed between the early histone gene loci on separate chromosomes and that the gene loci were directed to locate to more interior positions. Furthermore, these interactions were associated with the active transcription of the early histone genes. Thus, such dynamic interchromosomal interactions may contribute to the efficient synthesis of early histone mRNA during the morula stage of sea urchin development.

Reactivation of larval keratin gene (krt62.L) in blastema epithelium during Xenopus froglet limb regeneration.

Limb regeneration is considered a form of limb redevelopment because of the molecular and morphological similarities. Forming a regeneration blastema is, in essence, creating a developing limb bud in an adult body. This reactivation of a developmental process in a mature body is worth studying. Xenopus laevis has a biphasic life cycle that involves distinct larval and adult stages. These distinct developmental stages are useful for investigating the reactivation of developmental processes in post-metamorphic frogs (froglets). In this study, we focused on the re-expression of a larval gene (krt62.L) during Xenopus froglet limb regeneration. Recently renamed krt62.L, this gene was known as the larval keratin (xlk) gene, which is specific to larval-tadpole stages. During limb regeneration in a froglet, krt62.L was re-expressed in a basal layer of blastema epithelium, where adult-specific keratin (Krt12.6.S) expression was also observable. Nerves produce important regulatory factors for amphibian limb regeneration, and also play a role in blastema formation and maintenance. The effect of nerve function on krt62.L expression could be seen in the maintenance of krt62.L expression, but not in its induction. When an epidermis-stripped limb bud was grafted in a froglet blastema, the grafted limb bud could reach the digit-forming stage. This suggests that krt62.L-positive froglet blastema epithelium is able to support the limb development process. These findings imply that the developmental process is locally reactivated in an postmetamorphic body during limb regeneration.

Developmental changes in drug-metabolizing enzyme expression during metamorphosis of Xenopus tropicalis.

A large number of chemicals are routinely detected in aquatic environments, and these chemicals may adversely affect aquatic organisms. Accurate risk assessment requires understanding drug-metabolizing systems in aquatic organisms because metabolism of these chemicals is a critical determinant of chemical bioaccumulation and related toxicity. In this study, we evaluated mRNA expression levels of nuclear receptors and drug-metabolizing enzymes as well as cytochrome P450 (CYP) activities in pro-metamorphic tadpoles, froglets, and adult frogs to determine how drug-metabolizing systems are altered at different life stages. We found that drug-metabolizing systems in tadpoles were entirely immature, and therefore, tadpoles appeared to be more susceptible to chemicals compared with metamorphosed frogs. On the other hand, cyp1a mRNA expression and CYP1A-like activity were higher in tadpoles. We found that thyroid hormone (TH), which increases during metamorphosis, induced CYP1A-like activity. Because endogenous TH concentration is significantly increased during metamorphosis, endogenous TH would induce CYP1A-like activity in tadpoles.

Clustered Xenopus keratin genes: A genomic, transcriptomic, and proteomic analysis.

Keratin genes belong to the intermediate filament superfamily and their expression is altered following morphological and physiological changes in vertebrate epithelial cells. Keratin genes are divided into two groups, type I and II, and are clustered on vertebrate genomes, including those of Xenopus species. Various keratin genes have been identified and characterized by their unique expression patterns throughout ontogeny in Xenopus laevis; however, compilation of previously reported and newly identified keratin genes in two Xenopus species is required for our further understanding of keratin gene evolution, not only in amphibians but also in all terrestrial vertebrates. In this study, 120 putative type I and II keratin genes in total were identified based on the genome data from two Xenopus species. We revealed that most of these genes are highly clustered on two homeologous chromosomes, XLA9_10 and XLA2 in X. laevis, and XTR10 and XTR2 in X. tropicalis, which are orthologous to those of human, showing conserved synteny among tetrapods. RNA-Seq data from various embryonic stages and adult tissues highlighted the unique expression profiles of orthologous and homeologous keratin genes in developmental stage- and tissue-specific manners. Moreover, we identified dozens of epidermal keratin proteins from the whole embryo, larval skin, tail, and adult skin using shotgun proteomics. In light of our results, we discuss the radiation, diversification, and unique expression of the clustered keratin genes, which are closely related to epidermal development and terrestrial adaptation during amphibian evolution, including Xenopus speciation.

Rapid and efficient analysis of gene function using CRISPR-Cas9 in Xenopus tropicalis founders.

Recent advances in genome editing using programmable nucleases, such as zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALENs) and the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas system, have facilitated reverse genetics in Xenopus tropicalis. To establish a practical workflow for analyzing genes of interest using CRISPR-Cas9, we examined various experimental procedures and conditions. We first compared the efficiency of gene disruption between Cas9 protein and mRNA injection by analyzing genotype and phenotype frequency, and toxicity. Injection of X. tropicalis embryos with Cas9 mRNA resulted in high gene-disrupting efficiency comparable with that produced by Cas9 protein injection. To exactly evaluate the somatic mutation rates of on-target sites, amplicon sequencing and restriction fragment length polymorphism analysis using a restriction enzyme or recombinant Cas9 were performed. Mutation rates of two target genes (slc45a2 and ltk) required for pigmentation were estimated to be over 90% by both methods in animals exhibiting severe phenotypes, suggesting that targeted somatic mutations were biallelically introduced in almost all somatic cells of founder animals. Using a heteroduplex mobility assay, we also showed that off-target mutations were induced at a low rate. Based on our results, we propose a CRISPR-Cas9-mediated gene disruption workflow for a rapid and efficient analysis of gene function using X. tropicalis founders.

Cilia play a role in breaking left-right symmetry of the sea urchin embryo.

Left-right asymmetry of bilaterian animals is established during early development. In mice, frogs and fishes, the ciliated left-right organizer plays an essential role in establishing bilateral asymmetry, and leftward flow of extracellular fluid generated by ciliary motion results in Nodal activity on the left side. However, H(+) /K(+) -ATPase activity is also involved in the determination of left-right asymmetry in a variety of animals, and it has been thought to be an ancestral mechanism in deuterostomes. In sea urchin, the determination of the left-right asymmetry based on H(+) /K(+) -ATPase activity was already clarified, but it remains to be uncovered whether ciliary motion is involved in the left-right asymmetry of the embryo. Here, we show evidence that ciliary motion is involved in the establishment of left-right asymmetry of sea urchin embryo. Furthermore, we show that the initial cilia generated on small micromeres during the early stage of embryogenesis may be involved in this process. These results suggest that the cilia-mediated mechanism for the determination of left-right asymmetry may be acquired at the base of the deuterostomes.